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ABSTRACT 

Context. The nature of maser emission means that the apparent angular size of an individual maser spot is determined by the local 
, amplification process as well as by the instrinsic size of the emitting cloud. Highly sensitive radio interferometry images made using 

O MERLIN spatially and spectrally resolve water maser clouds around evolved stars. 
Aims. We used measurements of the cloud properties, around the red supergiant S Per and the AGB stars IK Tau, RT Vir, U Her and 
U Ori, to test maser beaming theory. In particular, spherical clouds are expected to produce an inverse relationship between maser 
O intensity and apparent size, which would not be seen from very elongated (cylindrical or slab-like) regions. 

Methods. We analysed the measured properties of the maser emission in order to estimate the saturation state. We analysed the 
variation of observed maser spot size with intensity and across the spectral line profiles. 

Results. Circumstellar masers emanate from discrete clouds from about one to 20 AU in diameter depending on the star. Most 
of the maser features have negative excitation temperatures close to zero and modest optical depths, showing that they are mainly 
unsaturated. Around S Per and (at most epochs) RT Vir and IK Tau, the maser component size shrinks with increasing intensity, 
although in some cases the slope is shallower than predicted, probably due to shape irregularities and the presence of velocity gradients 
within clouds. In contrast, the masers around U Ori and U Her tend to increase in size, with a larger scatter. 

Conclusions. The water masers from S Per, RT Vir and IK Tau are mainly beamed into spots with an observed angular size much 
smaller than the emitting clouds. The brighter spots at the line peaks are smaller than those in the wings. This suggests that the masers 
are amplification-bounded, emanating from clouds which can be approximated as spheres. Many of the masers around U Her and U 
Ori have apparent sizes which are more similar to the emitting clouds and have less or no dependence on intensity, which suggests that 
these masers are matter-bounded. This is consistent with an origin in flattened clouds and these two stars have shown other behaviour 
indicating the presence of shocks which could produce this effect. 

Key words. Masers - Stars: AGB - Stars: supergiants - circumstellar matter 

1. Introduction Unfortunately, maser geometry cannot be determined di- 
rectly. As noted in EHM92, from any given direction a spherical 

Maser regions are comprised of many bright, compact spots, maser appearfj identical to a cylindrical maser aligned along the 

£ . each with its own well-defined Doppler velocity. Geometry is Une of sight while the length of thifj equivalent cylinder is the 

Cd ; an important factor in controlling the brightness of these in- sphere , s diametel - its cross sect ion, the S p her ical maser observed 

dividual features. Elitzur, Hollenbach, & McKee (1992), here- ^ is controlled by the maser beaming and determined by the 

after EHM92, suggest that the main difference between the high amplification optical dept h ("amplification-bounded" maser). 

brightness temperatures of H 2 masers in star-forming regions The observed area of an amplification-bounded maser varies 

and their lower brightness counterparts in evolved stars is the inversely with the amplification along the line of sight— the 

geometry— the pumping scheme and physical conditions are stronger is the amplificatiorii the sma ii er is the observed area, 

similar in both but the former are characterised by planar shock Because the amplification decreases with frequency shift away 

geometry while the latter are three-dimensional structures that from Hne centre> the appearance of a sph erical maser varies 

can be modelled as spheres. The spherical maser radius is al- across the Une file . the observed area is smallest at line cen . 

ways limited because loss lines of the pump cycle become opti- ^ increasing toward the line wings (Elitzur 1990) In contrasti 

cally thick for sufficiently large dimensions, and the maser ther- the observed ^ of a cylindrical maser is its actual cross section 

mahzes. In contrast, there are no similar limitations on the size ( " matter _ bounded " maser ), independent of frequency. Therefore, 

of either disk or filamentary masers. In these cases the pumping the mafjer geometry can be determined indirectly by examin- 

cycle thermal photons can escape along a short dimension per- ing its observed size at different frequencies across the line pro- 

pendicular to the line of sight while the length along the latter file The cartoon Rg { illustrates these properties. The spherical 

can increase arbitrarily without any effect on the basic inversion douds on the lower left giye rise tQ the arri pii nC ation-bounded 

efficiency, allowing extremely high brightness temperatures. spatial profile shown at the upper left Some of the douds shown 
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at lower right are flattened and, if observed edge-on, matter- 
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Fig. 1. The lower part of this cartoon shows, on the left, spherical water 
vapour clouds in a CSE. On the right some clouds are flattened. The 
upper diagrams show the spatial profiles resulting from maser propa- 
gation through a sphere (left) and through a disc seen edge-on (right) 
The scale in terms of stellar radii is notional. The similarities with the 
named sources are explained later in this paper. 



bounded masing occurs, giving rise to the spectral profile at up- 
per right. 

In accordance with the EHM92 proposal, Genzel et al. 
(1981) report that observed sizes of individual H2O maser spots 
in W51M are roughly constant across the line profiles, and sim- 
ilar results were obtained for a number of other sources by R. 
C. Walker & J. M. Moran (private communication); the fre- 
quency independence of the observed sizes of bright masers in 
star-forming regions conforms with the matter-bounded geome- 
try expected behind shock fronts. 

On the other hand, stellar wind masers are expected to be lo- 
cal regions carved out of the outflow by the requirements of ve- 
locity coherence along the line of sight. The resulting geometry 
is in general fully three-dimensional and amplification-bounded, 
thus the observed area should increase away from the line cen- 
tre. Here we test this hypothesis for H2O masers in circumstellar 
envelopes (CSE) of evolved stars. 

H2O masers at 22 GHz are found around post-main- 
sequence oxygen-rich asymptotic giant branch (AGB) and red 
supergiant (RSG) stars, with surface temperatures of 2500- 
3500 K. AGB stars supporting water masers have radii of ~ 
0.5 - 3 AU and undergo mass loss at 10~ 6 - 10 5 M yr" 1 ; RSG 
are around an order of magnitude larger and more prolific. The 
masers occur within five to a hundred stellar radii (/?*) of the star. 
Individual features typically span 1-2 km s _1 and the observed 
size in a single 0. 1 km s _1 channel is typically a few milli-arcsec 
(mas). 




10 20 30 

Component size s (mas) 

Fig. 2. Histogram of the sizes of H2O maser components exceeding the 
observational error (s v > cr s ). The ordinate shows the number of com- 
ponents per size bin. The thickest black line shows the distribution for 
the epoch with the brightest maximum emission (/ milx , given in Table 2) 
for each star, with thinner and paler/bluer lines representing epochs with 
progressively fainter maxima. A few U Her and U Ori components with 
sizes > 40 mas are not shown. 



MERLIN 1 can produce detailed images of H2O maser shells 
around AGB stars within a few hundred pc, or RSG within a few 
kpc from Earth, without losing significant emission due to miss- 
ing interferometer spacings. We measured the sizes of individual 
maser spots by fitting Gaussian components, allowing us to in- 
vestigate the variations of maser size with intensity, position in 
the line profiles and other factors. The sample and observational 
results are described in Section 2. We summarise the evidence 
for maser clumping into discrete clouds and their properties and 
saturation state in Section 3. The variations of the observed 
size of maser components with intensity and with position in 
the spectral line profile of clouds are analysed in Section 4. 
This provides diagnostics for the two different maser geome- 
tries, explained in Section 5. Section 6 summarises compelling 
evidence that the variations of observed maser sizes across line 



1 The UK radio interferometer, operated by the University of 
Manchester on behalf of STFC 
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Table 1. Properties of the sample stars. The stellar velocity V* is given 
in the Local Standard of Rest (LSR) convention. The distances D used 
in our calculations are 2.3 kpc for S Per, 266 pc for the Miras, IK Tau, 
U Ori and U Her and 133 pc for the SRb RT Vir, in order to remain 
consistent with RYC99 and B+03. The values of R+ for the AGB stars 
were measured using IR interferometry at H and K bands. The radius of 
S Per is the average of values deduced from spectral fits. References in 
the table are A97 Alvarez & Mennessier (1997); C91 Chapman et al. 
(1991); C94 Chapman et al. (1994);C07 Cenarro et al. (2007); D87 
Diamond et al. (1987); H74 Humphreys (1974); K87 Kirrane (1987); 
L05 Levesque et al. (2005); M04 Monnier et al. (2004); M08 Mayne & 
Naylor (2008); N86 Nyman et al. (1986); 098 Olofsson et al. (1998); 
R06 Ragland et al. (2006); V07 Vlemmings & van Langevelde (2007); 
vL07 van Leeuwen (2007). 



profiles can be used to distinguish between matter-bounded and 
amplification-bounded masers. 

2. Observational measurements 

2.1. Observations 

We have carried out a program of MERLIN observations of 22- 
GHz H2O masers around evolved stars, from 1994 onwards. 
Four AGB stars (U Her, U Ori, IK Tau and RT Vir) and one RSG 
(S Per) have good enough visibility plane coverage to resolve 
the individual maser components. Their relevant properties are 
summarised in Table 1 . 

Observational parameters are given in Table 2. The cr rms 
noise varies greatly not only due to the different durations of 
observations but to the elevation, the peak (determining the ef- 
fectiveness of self calibration) and most of all the weather. The 
1994 observations and image analysis methods are described in 
detail by Bains et al. (2003) (B+03) for the AGB stars IK Tau, 
RT Vir, U Her and U Ori and in Richards et al. (1999a) (RYC99) 
for the RSG S Per. Similar procedures were used at later epochs, 
except that phase-reference sources were included after 1999 in 
order to improve position accuracy. The 1994 observations were 
made with channels of 0.105 km s _1 in a total span of 25 km 
s _1 , as was used for all epochs of U Ori, U Her and RT Vir. The 
spectra of IK Tau and S Per were found to have a larger velocity 
extent, so double the span was used for these two sources in sub- 
sequent epochs, necessitating a doubling of the channel width 
to 0.21 km s _1 . The image weightings and restoring beams were 
adopted to optimise the accuracy of component fitting, which de- 
pends mainly on the signal-to-noise ratio (RYC99) as long as the 
emission has a naturally Gaussian distribution. The masers are 
well-resolved by MERLIN; the total extent of the 22-GHz emis- 
sion regions is over a hundred mas, at a resolution of 10-20 mas, 
and many spectral channels contain multiple distinct patches of 
maser emission. 



RYC99 showed that we detect at least 85% of the S 
Per 22-GHz maser emission. Single-dish monitoring with the 
Pushchino telescope is available within a few months of some 
epochs for some of our sources. Lekht et al. (2005) detected 
slightly less flux than we did for S Per in 1994, and more in 1999. 
Rudnitskij et al. (2000) show that we detected a similar amount 
of maser emission from U Ori as they did in 1994 and 1999. The 
single-dish flux of RT Vir within a few months of our observa- 
tions (Mendoza-Torres et al. 1997; Lekht et al. 1999) was lower 
than the closest interferometry measurements, but both sets of 
observations show that it is highly variable, changing flux den- 
sity by a factor of 2 in a few months. 

2.2. Analysis of maser images 

Detailed results from the 1994 observations are given in RYC99 
and B+03 and the large-scale behaviour revealed by more recent 
observations will be analysed in a future paper. Here, we concen- 
trate on the properties relevant to understanding maser beaming, 
which are summarised in Table 2. 

Identifying individual velocity-coherent maser features in an 
expanding, turbulent outflow is a difficult task. To this end we 
developed a two-stage process to disentangle the maser struc- 
ture properly in velocity and space, first by characterising the 
spatial distribution of the emission channel by channel and then 
by examining the spatially-resolved spectral profiles. We start by 
introducing the first concept, of a maser component: Emission 
measured in a particular channel centred on velocity v, in a con- 
tiguous area on the plane of the sky centred on a given direction. 
The emission pattern of a 3D unsaturated maser in a given fre- 
quency interval has the Gaussian shape e~ (e/6v) , where 6 is dis- 
placement from the direction of the longest chord along the line- 
of-sight and 8 V is the beaming angle; the observed area is then 
A = n6l (Elitzur 1992). This property of maser emission allows 
us to obtain accurate measurements by fitting a two-dimensional 
Gaussian component to each patch of maser emission, in each 
channel, and deconvolving the restoring beam from the map area 
as described in RYC99. We selected components brighter than 4 
or 5<r rms noise (see RYC99 and B+03 for more details of re- 
jection of artefacts). The resulting fit determines the intensity /, 
position (x, y) and angular FWHM (full width at half maximum) 
s v (= 2 Vln2 6 V ) of each maser component. 

The uncertainty cr 5 in the component size is the sum in 
quadrature of the position uncertainties in the x and y directions. 
There is a lower limit of 0.1 mas to the resolvable component 
size (Richards 1997) due to dynamic range and calibration lim- 
itations. Between 1% and 40% of component sizes were upper 
limits and if cr, > j„,we took s v = cr s . Figure 2 shows the distri- 
bution of maser component sizes (restricted to s v > cr s ) for each 
source and epoch. The brightness temperature of a component is 



where / is in Jy, the wavelength A = 0.013 m, k B is Boltzmann's 
constant and the area A = (n/4- In 2)s 2 v is in steradian. 

Figs. 3-7 show the variations of s v across the line profiles at 
each epoch. Figs. 8, 12 and 15 illustrate the component distribu- 
tions in the most well-filled, symmetric CSE (S Per), the most 
sparse (U Ori) and a well-filled CSE with velocity asymmetry 
(RT Vir); see plot captions for details. These Figures show that 
maser components lie in series like pearls on a necklace, often 
curved or twisted. 
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Star Date Dur Bm cr rms r, r a Vi v a e NC NF / max / N Rt AVi /2 

(yymmdd) (hr) (mas) (mjyr 1 ) (AU) (km s' 1 ) (Jy) (AU) (km s' 1 ) 

SPer 1 940324 143 10 17 55 175 9l) 16 05l 1040 93 72 18^9 50 0.77 ± 0.29 

SPer 990110 13.0 10 8 39 175 8.0 22 0.56 689 100 47 12 ±6 56 0.91 ±0.34 
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UOri 
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20 


20 
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2.5 


5.5 


0.72 


94 


13 


44 


4.8 ±2.6 
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0.51 ±0.16 


UOri 


000410 
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18 


25 


7 


36 


2.0 


6.0 


0.65 


142 


31 
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2.7 ± 2.2 
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0.46 ± 0.23 


UOri 
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18 
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0.68 
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4.9 ±4.4 
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0.74 ± 0.23 
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13.7 
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14 
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9.5 
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2.3 ± 1.7 


18 


0.68 ±0.35 
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0.72 


212 
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3.9 ±3.5 
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0.62 ± 0.32 
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940415 
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10 


16 
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5.0 


16 
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25 
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11 
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41 
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1.0 ±0.8 


20 


1.05 ±0.42 


RTVir 


960524 


7.7 


12 


20 
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19 


3.5 


11 


0.81 


401 


51 
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0.9 ±0.7 


19 


0.88 ±0.31 


RTVir 


960612 


8.4 


12 


35 


5 


19 


3.5 


11 


0.81 


246 


38 


791 


0.9 ± 0.6 


19 


0.80 ± 0.25 



Table 2. The first 5 columns give the target, date, duration, resolution and sensitivity of each observation, r\, r„, Vi and v a are the inner and outer 
limits of the H2O maser shell, e is the logarithmic velocity gradient (Section 3.1.4). The total number of individual maser components (NC) and 
features (NF) made up of contiguous series of components are given, along with the maximum total intensity / max of the brightest spatially distinct 
component. / is the mean size of each feature, N m is the number for which a Gaussian spectral profile could be fitted at > 3cr and AV [/2 is the mean 
FWHM of these. All observations used 0.105 km s" 1 channel width apart from S Per epoch 9901 10 and IK Tau epochs 000520 and 010427, when 
0.21 km s -1 was used. 

l NC, NF and V are probably underestimates since the observational bandwidth was less than the likely velocity extent of emission. 



We introduce the second concept, of a maser feature: A se- 
ries of adjacent components within a beam width, spanning at 
least 3 velocity channels. Each feature was defined as possess- 
ing a single-peaked spectral profile, unless the emission was too 
faint to show a clear peak or the profile was truncated by blend- 
ing. All components were allocated to features consisting of col- 
lections of neighbours within a beam width in successive chan- 
nels. In cases of ambiguity, component allocation was guided 
by following the smoothest increments of position and flux den- 
sity from channel to channel. A few, isolated components were 



rejected as artefacts, being very faint and/or in the position of 
beam side lobes. 

An example feature is enlarged in Fig. 9 (region E in Fig. 8). 
This shows that the smooth line profile corresponds to a curved 
position-velocity gradient. The angular FWHM d of the feature 
is given by the separation of components with intensities closest 
to half maximum intensity, marked by the short fat red line. The 
largest angular size I (uncertainty erf) between its most widely 
separated components is shown by the long red line in Fig. 9. 
The interpretation of d is discussed further in Section 3.2.4. 
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Fig. 4. The velocity profile of U Ori, see Fig. 3 for details. 
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Fig. 5. The velocity profile of U Her, see Fig. 3 for details. 



We consider the features thus identified as the fundamental 
building blocks of the maser source. We fitted the spectral profile 
of each feature, i.e. the component intensity I(v) as a function of 
v, with a Gaussian curve. This is an approximation to the line 
shape of unsaturated masers, adequate to half maximum. This 
fitting yields the line FWHM AV1/2, and the corresponding un- 
certainty cr A y 1/2 , as well as the central velocity, the peak model 
intensity and their uncertainties. The relevant measurements are 
given in Table 2, along with the number of features at each epoch 
with fits better than 3a-. 
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Fig. 6. The velocity profile of IK Tau, see Fig. 3 for details. 



We show additional enlargements of features in Figs. 10, 11, 
13, 14, 16 and 17 as clear (but otherwise unexceptional) illustra- 
tions of different relationships between s v and intensity or posi- 
tion in the line profile, discussed in Section 4. Figs. 10 and 11 
show multiple spectral peaks in spatially adjacent features which 
probably emanate from a single physical region, discussed in 
Section 3.1.1. 

In summary, the individual maser components represent sep- 
arate patches of emission as sampled in the observed velocity 
channels. The components are found to be clustered in position 
and velocity, often forming smooth Gaussian spectral profiles, 
termed features. Individual features thus defined are the objects 
to be compared with maser theory. Table 2 gives the number of 
components and features in each source. 

3. Properties of maser clouds 

3. 1 . Dense clumps in the stellar winds 

The evidence for the origins of 22-GHz H2O maser emission in 
discrete, density-bounded clouds, as part of a fragmented, ex- 
panding wind, is summarised in this section. RYC99 and B+03 
explained in more detail why we think that each feature (or occa- 
sionally a series of a few features, in the RSG) corresponds to a 
physically discrete, density bounded cloud, not solely a 'spooky' 
accident of the velocity field or turbulence (Strelnitski 2007). 

3.1 .1 . Maser components form distinct features with 
Gaussian spectral profiles 

The relationship between component positions and velocities is 
illustrated for S Per (Figs. 8 and 10), U Ori (Figs. 12 and 13) and 
RT Vir (Figs. 15 and 16). Many features show an orderly (but not 
necessarily linear) distribution of components, with an internal 
gradient of velocity with position. These gradients do not show 
any large-scale alignment. The average line widths AV1/2 for the 
feature spectral profiles for each star are given in Table 2. The 
gas temperature T in the CSE of AGB stars is likely to drop from 
around 1000 K at r\ to 400 K at r Q , depending on the thickness 
of the shell (Zubko & Elitzur 2000). 
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Fig. 7. The velocity profile of RT Vir, see Fig. 3 for details. 96a to f 
indicate the six successive epochs of observation in 1996. 



The temperature at distance r from the star can be approx- 
imated as T oc r~° , so the thermal Doppler line width AVd oc 
Vr oc r~°- 2 . AVd lies in the range 1.6-1.2 km s _1 between r\ and 
r . The weak temperature dependence means that even if RSGs 
and their CSEs are slightly hotter than AGB stars (Table 1), 
this has little effect on maser amplification paths. AVd exceeds 
AV1/2, so maser amplification is likely to take place through the 
entire depth of each feature. 

Features around AGB stars are almost always distinct in po- 
sition or velocity (or both) and we conclude that they correspond 
to separate clouds. A few features around RSGs form groups 
such as the example in Fig. 8 region A. This appears to be a 
single cloud containing 5-6 multiple, adjacent spectral peaks 
(Fig. 11), i.e. multiple features in our terminology. These large 



E -20- 




80 60 40 20 -20 -40 -60 -80 
Right Ascension offset (mas) 

Fig. 8. Each symbol represents an H 2 maser component in a single 
channel observed towards S Per in 1994 (o) or 1999 (A). The diameter 
is proportional to the square root of total intensity. The colour is propor- 
tional to the velocity as shown. The red boxes enclose features which 
are shown in more detail in Figs. 9, 10 and 11. 



clouds have a velocity span > AVd and RYC99 showed that their 
velocity gradients, rather than their total depths, probably limits 
their maser amplification. The majority (80% - 90%) of features 
around RSGs appear to emanate from separate clouds. 

3.1 .2. The angular extent of features gives the clouds 
physical sizes 

Emission from successive velocity channels sampling a quies- 
cent sphere would be centred on a common point on the sky. 
Alternatively, a long thin line of observed maser components 
would be produced by a spherical cloud with a linear internal 
velocity gradient. The patterns we observe could be produced by 
spherical clouds with curved velocity gradients. It seems likely 
that MERLIN is detecting all the maser emission (Section 2.1) 
so / is close to the full extent of the maser features with signifi- 
cant position-velocity gradients. In other cases, I is a lower limit, 
but individual clouds cannot exceed about 20 / in any direction 
or they would not fit inside the observed maser shell inner and 
outer radii. 

The clouds are not necessarily spherical but could be amoe- 
boid or have any other shape. Unfortunately we cannot measure 
all three dimensions of each cloud directly so we use I (con- 
verted to units of length at the appropriate distance for each star) 
as the physical diameter of a feature in any direction, indicating 
the most likely biases this causes where relevant. The analysis in 
this paper (e.g. Section 5) will provide further clues to the shapes 
of clouds. 



3.1 .3. Maser clouds persist for many months or years 

The distribution of components within some features forms dis- 
tinctive patterns which can be recognised (with some distortion) 
at multiple epochs. Close similarities can be seen in some in- 
stances, such as the RT Vir feature shown in Fig. 16 (one of 1 1 
features matched at all 6 epochs of observation taken over 10 
weeks). The S Per features seen in Fig. 8 region A, also shown 
in Fig. 10, appear to have rotated about 60° anticlockwise during 
5 years. The epochs were aligned using the centres of expansion 
(found as described in RYC99), combined, relative uncertainties 
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Fig. 9. Properties of the maser feature E in Fig. 8, with the same col- 
oring scheme. Top: Feature line profile. The symbol diameters are pro- 
portional to s v (two are very small); the error bars show the size un- 
certainty cr s . Bottom: Detailed map of the components comprising the 
feature. The diameter of each symbol is 10% of s v and the error bars 
represent the position offset uncertainties. The thin black line joins the 
components in velocity order. The long thin red line denotes the largest 
angular size / of the feature. The short fat red line shows the separation 
of the components with flux densities closest to the half-maximum of 
the feature, d. The expanded-scale inset (total size 1.7 x 1.2 mas, ticks 
at 0.2 mas intervals) shows the area in the pink box, which contains 
the components brighter than half maximum, with component size er- 
ror bars. The brightest component is outlined in red, also marked on the 
Top plot. 



~ 6.5 mas. Proper motion studies of the AGB stars show that 
masing from their H2O clouds can be tracked for between a few 
months and a couple of years (Richards et al. 1999b; Yates & 
Cohen 1994). The larger RSG clouds can be tracked for 5 to 
10 or more years (Richards et al. 1996; Richards et al. 1998; 
Murakawa et al. 2003). 

3.1 .4. Water maser clouds are much denser than their 
surroundings 

The H2O masers propagate in approximately spherical shells 
with inner and outer radii of r\ and r Q , across which the expan- 
sion velocity increases from v\ to v . Multi-epoch and proper 
motion measurements suggest that the H2O maser shells are 
approximately spherically symmetric (Murakawa et al. 2003, 
Richards et al. 1999b). This allowed Murakawa et al. (2003) 
(their equation 1) to solve the three-dimensional structure of the 
circumstellar envelopes and find the distance of each cloud from 
the star, r. 

The wind acceleration through the shell is parametrised by 
the logarithmic velocity gradient e = log(ii /uj)/ log(r /n). 
Alcock & Ross (1986) proposed that if initially spherical, dusty 
clouds are driven away from a star by radiation pressure on 
grains, at a constant velocity, their diameter in the radial, out- 
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Fig. 10. The angular size, distribution and position uncertainties of the 
maser components making up the boxed features A-D around S Per in 
Fig. 8, as measured in 1994 (upper and middle plots, circles) and the 
features in A as measured in 1999 (lower plot, triangles). The same 
velocity colour scheme is used. Region A contains 6 features in 1994 
and 5 in 1999; the other regions contain 2 features, as can be seen from 
the spectral peaks in Fig. 11. The diameter of each symbol is 10% of 
s v . The black lines join the components in velocity order. The position 
uncertainties are shown. 



flow direction will be unchanged but their diameter in the tan- 
gential direction will increase in proportion to distance from the 
star. If the wind is accelerating, then the tangential and radial 
axes (with respect to the star) of a cloud expanding under radi- 
ation pressure will grow as r/r; and (r/r;) 6 , respectively. If the 
clouds are spherical at r u their tangential/radial aspect ratio will 
evolve as (r/>i) 1_e . The cloud number density, relative to that at 
the inner radius, is given by n{r) = n r . x (r/r;)~ (2+e) . 

The inner radius r\ is determined by the distance from the 
star where the wind density falls below the collisional quench- 
ing density for the 22-GHz H2O maser (Cooke & Elitzur 1985; 
Cohen 1987), at a number density of n(r\) « 5 x 10 15 m -3 
(a: 5 x 10 9 cm -3 ). The mass loss rate M which would be needed 
to produce n(r;) is one or two orders of magnitude greater than 
the mass loss rates measured from IR and CO observations (see 
references in RYC99 and B+03). 

This apparent contradiction is overcome by considering the 
volume filling factors of the H2O maser clouds. The number of 
clouds NC and their sizes I are given in Table 2 and the filling 
factors are given by A^C x / 3 /(r 3 - r 3 ). S Per has the highest fac- 
tor of 1.5% (RYC99), RT Vir and IK Tau have factors of about 
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Fig. 11. Top to bottom: The line profiles of features seen around S Per 
in 1994 (circles) in the regions labelled A, B, C and D in Fig. 8, and 
in 1999 (triangles) in the region labelled A in Fig. 8. The diameters of 
the circles/triangles and of the crosses are proportional to the observed 
angular FWHM s v of the maser components and to <x„ respectively. 



0.4% and U Ori and U Her have factors of about 0.1% (B+03). 
This leads to realistic values of M if the clouds are surrounded 
by gas at ~ 1/50 of the cloud number density. This implies that 
the actual average feature volume per CSE must be no more than 
a factor of about two greater than that based on the average mea- 
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Fig. 12. Each symbol represents an HiO maser component in a single 
channel observed towards U Ori, at the 4 epochs indicated. The diam- 
eter is proportional to the square root of total intensity. The colour is 
proportional to the velocity as shown. The numbered features are shown 
in more detail in Figs. 13 and 14. 
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Fig. 13. The angular size and distribution of the maser components in 
the U Ori features marked on Fig. 12, with the same colouring scheme. 
Each panel shows a single feature. The diameter of each symbol is 10% 
of s v and the error bars represent the position offset uncertainties. The 
line joining the components shows the velocity order. 



sured /, or the mass in features would exceed the total mass loss 
rate. 



3.2. Maser optical depth and saturation state 

In this Section, we use our measurements to provide estimates 
of the parameters determining the maser optical depths, which 
leads us to infer that the masers are mostly unsaturated. We pro- 
vide estimates of the uncertainties in our measurements and the 
underlying assumptions; although these are large, different ap- 
proaches give consistent results. 



3.2.1. Estimating the population inversion 

An = «2 - "l is the difference between the number densities 
per substate of the upper and lower levels involved in the maser 
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eters of the circles and crosses are proportional to the observed angular 
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transition, 6i6 and 523 of ortho-water (de Jong 1973). Maser am- 
plification requires An > 0. 

We start by assuming that the emission rate is approaching 
the pump rate but not limited by it, following the assumptions 
outlined in RYC99 (section 4.5 and references therein), which 
explains why this is reasonable, at least for the brighter masers. 
The population difference is given by equation 6.32 of Reid & 
Moran (1988). RYC99 (their equation 16) derived an approxi- 
mation in terms of measurable quantities; omitting constants this 
is: 



An oc 



AVi /2 rB 

YA 2 l 



(2) 



for spherical clouds. Y is the decay rate from the excited state, 
initially taken as 1 s _1 (Reid & Moran 1988). 

We obtained mean values of 10 3 < An < 4 x 10 4 irT 3 
(0.001 < An < 0.04 cm 4 ) for S Per, U Ori, UHer and IK 
Tau. RT Vir had consistently high values, 10 4 < An < 9 x 10 s 
irT 3 (0.01 < An~ < 0.9 ctrT 3 ). 

The population of the lower state n\ was estimated using 
(«i/wh 2 o)("Hio/"h 2 )«h 2 - We took the fractional number density 
of H2O with respect to H2 («h 2 o/"h 2 ) as 2x 10~ 4 , where the frac- 
tion of H2O molecules in the n\ state is 5 x 10~ 3 . The actual val- 
ues are likely to differ by a factor of 2 or more, both between and 
within sources. We take «h 2 as 65% of the total number density 
n, and estimate n as a function of r as outlined in Section 3.1.4. 

If the pump rate is much higher than the stimulated emission 
rate, An is a lower limit. If the masers are significantly saturated 
then the present assumptions about the geometry are affected and 
An could be overestimated. Errors in r are unlikely to exceed 
30%, producing less than a factor of ~2 uncertainty in n, and 
hence n\. An will be overestimated if Y or the cloud depth are 
underestimated. The average cloud depth is unlikely to be more 
than about twice I and no individual cloud depth can exceed ~ 
20/, as explained in Section 3.1.4. Y is unlikely to be that much 
greater than our assumed value, so An is very unlikely to be more 
than a factor of 20 too large. 

The mean percentage population inversion (weighted by 
measurement errors), 100 X Anj{n\ + 112), was ^ 1% for all 



sources apart from five RT Vir epochs and one IK Tau epoch 
which reached between 1% and ~ 10%. Selecting brighter fea- 
tures, with the brightness temperature at line peak r B o > 10 9 K, 
gave a higher mean population inversion of a few percent for S 
Per, U Ori and U Her at most epochs with meaningful data and 
a few tens of percent for one IK Tau epoch and RT Vir at most 
epochs. 

3.2.2. Maser excitation temperature 

The maser excitation temperature is defined as 



T x = (hv/k B )/ln[ni/(An + m)] 



(3) 



where h and k B are Planck's and Boltzmann's constants respec- 
tively. For v = 22.235 GHz, hv/ks ~ 1 K. We derived T x at each 
feature peak using our values of An and n\. The error- weighted 
mean T x for each epoch (given in Table 3) was -5 to -0.2 K. S 
Per, U Her and RT Vir have T x ^ - 1 K, as expected for the con- 
sistently brighter sources. The logarithmic relationship means 
that even if An is decreased by a factor of 20, most of the values 
of T x are decreased to not less than -20 K, apart from for U Ori 
where it approaches -100 K, although based on few significant 
measurements. 

3.2.3. Maser optical depth 

The maser optical depth for the brightest part of each cloud is 
related to the brightness temperature by 



r = \n(T BQ /\T x \). 



(4) 



The mean (T^o) and extrema of Tbo for each source are given 
in Table 3. The maxima may be underestimated, if the brightest 
peak component size is less than the minimum resolvable size. 
The numbers of features, NF, detected at each epoch (Table 2), 
do not have any significant relationship with 3cr rms (except pos- 
sibly for IK Tau), suggesting that the minima are not sensitivity- 
limited. 

We derived t using our estimated values of T x . Table 3 gives 
the mean values f for each epoch. The dispersion includes the 
scatter of measured values but not the uncertainties in the un- 
derlying assumptions. 7 lies between 8 and 16 for S Per, U Her, 
IK Tau and the epochs of U Ori with meaningful data. RT Vir 
has consistently higher values reaching f — 25. If An is overes- 
timated by a factor of 20, r would be reduced by ^ 5. If An is a 
lower limit, T x and t are also lower limits. 

3.2.4. Maser beaming angle 

The full angular size of a cloud and the apparent angular FWHM 
of maser features are approximated by / and d (as marked in 
Fig. 9). The small size of d means that the relative measure- 
ment errors are large and its interpretation depends on the in- 
ternal cloud structure. Figs 9, 10 and 16 show that the smallest 
components sampling the velocity channels near the line peak 
of bright features are often more closely-spaced on the sky than 
the fainter components in the line wings. If the maser features 
came from completely quiescent spheres, d would be vanish- 
ingly small. We calculated the average ratio d/so for each object 
at each epoch, where so is the observed size of the single maser 
component closest to the peak of the line profile fitted to each 
feature. This ratio is in the range 0.3 to 2.4, error-weighted mean 
0.9 + 0.2. This is consistent with a systematic velocity gradient 
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within features of similar magnitude to the thermal velocity dis- 
persion. 

Subject to these caveats, we can obtain order-of-magnitude 
estimates of the beaming angle O e st = d 2 /l 2 . RYC99 obtained 
an average of 1.5 + 0.8 x 10~ 3 sr for the 5 brightest features in S 
Per, which suggested that they were on the verge of saturation. 
We obtained error-weighted mean values per epoch, £2 est , in the 
ranges (0.8 - 6) x 10~ 3 sr for S Per, (0.6 - 2) x 10 3 sr for U Her 
and (0.5 - 10) x 10 3 sr for RT Vir. U Ori and IK Tau gave larger 
values, of (3 - 100) x 10~ 3 sr and (10-50) X 10~ 3 sr, respectively. 

Vlemmings et al. (2005), following Nedoluha & Watson 
(1992), show that masers remain unsaturated as long as T^O. < 
10 10 K sr, for F — 1 s _1 . We obtain average values of Tsoilest < 
10 9 K sr for all sources apart from RT Vir, where 7Bo£2est ~ 10 11 
K sr. However, the rapid variability of RT Vir (Richards et al. 
1999b; Lekht et al. 1999) suggests that its masers are also largely 
unsaturated, possible if F > I. Vlemmings & van Langevelde 
(2005) estimate T B Q ~ 10 10 - 10 11 for S Per from line pro- 
file analysis and deduce that the masers are mostly unsaturated, 
although their VLBA observations only detected the brighter, 
more compact emission. 

3.2.5. Maser line narrowing 



The maser FWHM linewidth is given by AVi/ 2 = A V D / V 1 + r 
(Elitzur 1992) for unsaturated amplification of ambient thermal 
radiation. This predicts 0.25 S5 AV\/2 ^ 0.4 for r in the range 
shown in Table 3 and AVd between 1.2-1.6 km s . The aver- 
age measured values of AV1/2, given in Table 2, are greater than 
predicted, whilst narrower than AVd- Approaching saturation 
causes the line profile to re-broaden. Alternatively, the observed 
linewidths may be explained by the velocity gradients within the 
clouds. Vlemmings & van Langevelde (2005) show that, in the 
presence of a velocity difference along the maser path of up to 
AVd, the unsaturated maser line width is up to about double the 
width at constant velocity. 



3.2.6. Predominantly unsaturated masers 

A large value of An provides stronger population inversion 
which provides conditions for bright, unsaturated emission (see 
Reid & Moran (1988), their fig. 6.2), and a negative T x of small 
magnitude, approaching zero from below. However, if the maser 
approaches saturation, An becomes small and T x approaches -00 
as the maser emission rate becomes close to the pump rate. The 
small-magnitude negative values we obtained for T x suggest that 
even the most intense H2O masers in these objects are unlikely 
to be strongly saturated. Our estimates of beaming angle and line 
narrowing also rule out strongly saturated water masers. 

We investigated the relationship between excitation temper- 
ature and peak brightness temperature for each source by eval- 
uating the slope A of the expression log|T x | oc AloglBo, for 
measurements more significant than than 3<x. This is an inverse 
relationship at every epoch. RT Vir has the shallowest slopes, av- 
eraging -0.6 ± 0.2; the average slopes for the other sources are 
U Her -0.8 ± 0.1, S Per -0.9 ± 0.1 and IK Tau -1.0 + 0.2. U Ori 
has the least significant data, giving -0.6 ± 0.5. This shows that 
the brightest masers are associated with the smallest-magnitude 
values of T x , consistent with unsaturated emission. 

The maser optical depth t is equivalent to the product of 
the unsaturated absorption coefficient kq and the path length /. 
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Fig. 15. Each symbol represents an H 2 maser component in a single 
channel observed towards RT Vir, epoch 960524. The diameter is pro- 
portional to the square root of total intensity. The colour is proportional 
to the velocity as shown. The boxed feature is shown in more detail in 
Figs. 16 and 17. 



We can deduce from the arguments of EHM92 2 that a spherical 
maser begins to saturate when r ^ 17 and becomes fully satu- 
rated when t ^ 100. The values of t which we obtain are well 
within the unsaturated or slightly saturated regimes. All these re- 
sults imply that circumstellar 22-GHz H2O masers are predomi- 
nantly unsaturated. 

4. The observed size of maser components 

The brightest components near the line centres have more com- 
pact sizes in many of our sources, with larger components in 
the line wings. Figures 3-7 show the variations of s v across the 
line profiles at each epoch. s v increases significantly from the 
line centres to the wings of almost all S Per features (apart from 
the faintest, with larger uncertainties). Most RT Vir features also 
show an inverse variation of s v with intensity, although some 
components near the peaks are anomalously swollen. In contrast, 
many U Ori and U Her features have no apparent correlation be- 
tween s v and intensity. IK Tau shows a mixture of behaviours. 
Figs. 9, 11, 14 and 17 show the line profiles of the selected fea- 
tures labelled in Figs 8, 12 and 15. These are examples of differ- 
ences in appearance which are characteristic of the three sources. 
U Her is similar to U Ori and IK Tau is somewhat similar to RT 
Vir. 

Figs. 10, 13 and 16 show s v as a function of position for the 
selected features (at 20% of the actual measured size of s v , in 
order to avoid excessive overlaps). The error bars represent the 
position errors. These plots show that s v /l ranges from «: 1% 
for the peak S Per and RT Vir components, to > 50% for the 
components close to their line edges, whilst U Ori shows a less 
ordered scatter. 

In the following sections, we compare the observed rela- 
tionships between the sizes of individual maser components and 
their intensities, with the models developed by EHM92. We ten- 
tatively investigate the effects of the sizes and velocity gradi- 



2 note that EHM92 use / as the radius of a cloud whereas we use it as 
the diameter 
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Fig. 16. The angular size and distribution of the maser components in 
the RT Vir feature boxed in Fig. 15, with the same colouring scheme. 
Each panel shows a single feature. The diameter of each symbol is 20% 
of s v and the error bars represent the position offset uncertainties. The 
green inset shows an enlargement of part of the 960524 feature. 



ents of the parent clouds. A fully quantitative analysis is difficult 
due to the difficulties in measuring the three-dimensional struc- 
ture of the clouds nor of the beaming angle. Nonetheless, we 
are able to distinguish between the predictions for amplification- 
and matter-bounded masers and identify which aspects of maser 
behaviour are well-described by current theories and which con- 
found them. 



4.1. Spherical, uniform clouds 

We start by comparing our results with the simplest model, 
for spherical, density bounded clouds. EHM92 show that 
these would produce amplification-bounded masers wherein the 
beaming angle should increase with frequency offset from the 
line centre. Initially, we assume that all the clouds are the same 
size and all are spherical (i.e., that / is constant within each 
source). 
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Fig. 17. The line profiles of RT Vir features marked on Fig. 15. The 
diameters of the circles and crosses are proportional to the observed 
angular FWHM s v of the maser components and to cr s , respectively. 



4.1 .1 . Component size variation with intensity: the model 

The observed size s v represents the angular FWHM of the distri- 
bution of the intensity of maser emission from water molecules 
in a particular cloud in the velocity interval sampled by a single 
channel. Using expressions from Elitzur (1992) for the effect of 
the maser amplification process in spherical clouds, the beaming 
solid angle is given by Q = jt/Q.51kq v , for unsaturated emission 
from a sphere of diameter /, where the unsaturated absorption 
coefficient for that velocity channel is kq v . (If the emission is 
significantly saturated the beaming angle is reduced by a factor 
of r s /0.5l where r s is the radius of the unsaturated core.) The 
observed area is given by A — {ji IA\n2)s 2 v = Q(0.50 2 , so the 
apparent (unsaturated) component size should be given by 



si oc 0.5/M),. 



(5) 



The intensity of unsaturated emission from a single veloc- 
ity channel though a cloud of depth I is given by / w Soe* * 1 , if 
/ » So, the input radiation (generally the ambient local thermal 
radiation at 22 GHz). This is clearly the case in CSEs where the 
background radiation is undetectable by interferometric obser- 
vations on the scale of the masers. The stellar continuum could 
be amplified by a cloud along the line of sight, but there is no 
evidence for exceptionally bright blue-shifted peaks in our data. 
The volume filling factor is S5 1 .5% so cloud overlap is not likely 
to be significant. The maser component size is related to the 
emission intensity by I /So « e ^ so (assuming In / » In So) 



(6) 



log(s) oc log(V0.50 + alogOn/) 



4.1 .2. Relationships between observed component size and 
intensity 

We test the simple, spherical model for amplification-bounded 
masers (from clouds of constant size) by plotting log s against 
log(ln/). Equation 6 predicts that this should give a slope of a = 
-0.5. Figs. 18 - 23 show least-squares error-weighted fits to the 
data for all sources. The values of a are given in Table 3. a is 
negative, showing that the apparent component size shrinks with 
increasing intensity, at both epochs for S Per, at 6/7 epochs for 
RT Vir, at 2/3 epochs for IK Tau and at 1/4 epochs for U Ori. 
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Table 3. Properties of H2O masers. All temperatures are in K. The brightness temperature Tbo was derived using Equation 1. Tbo is the mean value 
for the peaks of all features in each source, given along with the minimum and maximum values for the feature peaks. The mean values of the 
maser excitation temperature T x and optical depth t were derived as described in Section 3.2, which explains the uncertainties and dispersions; the 
given values are probably lower limits, a parametrises the relationship between component size and intensity for spherical, same-size clouds, as 
defined in Equation 6. N/? is the number of components for which the line profile fitting was better than 3<r, s v > 3tr s , s v > sq and the component 
lies within a feature velocity span of < AV D . f3 parametrises the relationship between component size, position in the line profile and the effect of 
the internal velocity gradient see Equation A.2. Numbers in italics indicate unreliable results where < 30% of the data for that epoch meet these 
requirements. 



Figs. 22 and 23 show that RT Vir most closely exhibits the 
behaviour predicted for spherical clouds, with a mean value of 
a of -0.5 for all 7 epochs. The behaviour is fairly consistent 
between epochs. 

S Per shows a < -0.5 for both epochs, as plotted in Fig. 18. 
The uncertainty in s v is inversely proportional to I and hence 
cr s is greater for fainter components, biasing the fit towards the 
contribution of the brighter components. These are most likely 
to be approaching saturation, causing a steeper shrinkage of O 
with increasing maser amplification. There are sufficient data for 
S Per to examine the properties for subsets. Selecting the fainter 
components only, with I < 0.6 Jy, gives a = -0.6 and -0.4 
for the 1994 and 1999 data respectively, close to the predicted 
values for unsaturated emission. This suggests that the steeper 
slopes measured for the whole data sets are due to saturation of 
the brighter components. This is consistent with the conclusion 
of RYC99 that the brightest features are approaching saturation. 

IK Tau also produces a < -0.5 in 2000 and 2001, but in- 
spection of Fig. 21 suggests a similar explanation as for S Per, 
although attempting to fit to sub-sets of the data produces a 
very large scatter. The apparent slope of a > in 1994 is due 
to a small group of bright components showing anomalous be- 
haviour; if these are excluded by selecting I < 7.5 Jy, a slope of 
-1.6 is obtained, whilst for I < 0.75 Jy, a = -0.6, but with low 
significance since there are few points with accurate measure- 
ments. 

U Ori and U Her have a > for all but one epoch and Figs 19 
and 20 show very large scatters. This shows that most of their 
maser regions do not behave as spherical clouds, discussed fur- 
ther in Section 5.2. 




ln(l) 



Fig. 18. S Per maser component FWHM s v as a function of the natural 
logarithm of the intensity / coded per epoch as in Fig. 19. 



4.2. Cloud size variations and velocity gradients 

The observed size of clouds within each CSE varies by a fac- 
tor of two or more (Table 2). We investigated whether this can 
be included meaningfully in our models and attempted to al- 
low for the effects of velocity gradients within clouds. This is 
described in Appendix A. In brief, we modified Equation 6 to 
include the change in path length through a cloud with velocity 
offset. Fitting this to the data gave slopes roughly consistent with 
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Fig. 19. U Ori maser component FWHM s v as a function of the natu- 
ral logarithm of the intensity /, for components where s v > 3cr s . The 
magenta (hollow), cyan (solid), orange (cross) and grey (triangle) sym- 
bols represent 1994, 1999, 2000 and 2001 results, respectively. The cor- 
responding error-weighted least-squares fits (to all data) are shown as 
solid, dashed, dot-dashed and dotted lines. 




Fig. 20. U Her maser component FWHM s v as a function of the natural 
logarithm of the intensity /, coded per epoch as in Fig. 19. 



the expected relationship between component size and position 
in the line profile, for those sources which showed the behaviour 
expected from spherical clouds (i.e. a < in Equation 6). 
However there were very large scatters, unsurprising because 
we are measuring velocity along the line of sight but angular 
separation perpendicular to this. In reality the clouds are not 
perfectly spherically symmetric and we have insufficient con- 
straints on the three dimensional structure of individual clouds 
to improve these results. Nonetheless, insofar as the results in 
Appendix A are meaningful, they are consistent with the results 
in Section 4.1.2. 

4.3. Confirmation and limitations of the predictions for 
amplification-bounded masing 

To sum up these results, the measured size of maser compo- 
nents has an inverse relationship with intensity (Section 4.1) 
for S Per, RT Vir and most epochs of IK Tau, as predicted for 
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Fig. 21. IK Tau maser component FWHM s v as a function of the natural 
logarithm of the intensity /, coded per epoch as in Fig. 19. 
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Fig. 22. RT Vir maser component FWHM s v as a function of the natural 
logarithm of the intensity /, observed in 1994. 




Fig. 23. RT Vir maser component FWHM s v as a function of the natural 
logarithm of the intensity /, observed in 1996. 
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amplification-bounded masers. In some of these cases, the com- 
ponent size shrinks more steeply than expected. This is partly be- 
cause results are biased towards the brightest components, which 
may be approaching saturation. The remaining sources, U Ori 
and U Her, show a mixture of behaviours and may be matter- 
bounded, see Section 5. This first, simple investigation, of the 
relationship between component size and intensity, assumed that 
all the features in a given source/epoch arise from identical, qui- 
escent spheres. 

We next attempted to investigate the behaviour of compo- 
nents within features, with respect to the component at the peak 
of each line, for masers meeting the amplification-bounded crite- 
rion of Sv > sq (Appendix A). We derived a semi-empirical rela- 
tionship which included the effects of velocity gradients. The av- 
erage ratio of the size of a component with respect to the size of 
the line peak component does increase with displacement from 
the line centre, but more weakly than predicted. 

EHM92 note that the beaming angle has a linear (rather 
than quadratic) dependence on the ratio between the observed 
size and the depth of an emission region, when these are deter- 
mined by orthogonal velocity gradients. However, further work 
is needed to explore application of this model quantitatively to 
our data due to the probable unsaturated state and irregular ve- 
locity gradients. In future, we hope to use proper motion studies 
which could establish the relationship between the velocity gra- 
dients parallel and perpendicular to the line of sight. e-MERLIN, 
will have at least double the sensitivity of MERLIN in the same 
width velocity channel, and allow sampling of any velocity span 
at any desirable resolution. 

5. Beaming and the properties of masing clouds 

We have compared maser properties with the predictions of 
maser beaming models. We now examine the implications for 
the geometry of the individual masing features and the condi- 
tions in the CSEs as a whole. Table 2 shows that e > 0.5, and 
r Q ~ (3 - 4)r;. Applying the expressions given in Section 3.1.4 
suggests that if clouds are initially spherical, expanding radially 
away from the central star under radiation pressure alone, their 
tangentiakradial aspect ratio will remain < 2. Some of our re- 
sults are consistent with this scenario (Section 5.1) but in other 
sources some of the clouds may be more flattened (Section 5.2). 

5.1. Amplification-bounded masers 

The measured size s v of H2O maser components decreases with 
increasing maser intensity, for emission from S Per, and, for 
the majority of epochs, from RT Vir and IK Tau (Section 4). 
The relationships roughly follow the expected relationship for 
maser beaming from unsaturated, spherical masers (Equation 6). 
EHM92 call such masers "amplification bounded" since the am- 
plification process controls the observed size. 

The beaming angle should increase with the frequency off- 
set from line centre within an amplification-bounded maser fea- 
ture. The relative observed component size, \n(s v /so), is ex- 
pected to be proportional to the velocity displacement from the 
line peak; we also attempt to allow for the observed velocity 
gradients using a semi-empirical relationship (Equation A.2). 
Qualitative agreements are obtained for sources/epochs where 
amplification-bounded masing dominates but the slopes are shal- 
lower than predicted, probably due to the over-simplified nature 
of the model and the non-linear velocity gradients. 

One aspect of the behaviour of RT Vir is anomalous. Figs. 7 
and 17 show that, in some of the brightest clouds, the compo- 



nents in the line wings appear relatively large, shrinking at first 
as expected towards the centre of the line, but then the observed 
sizes re-broaden near the peaks of the spectral features. This 
behaviour is very variable, e.g. the smallest components of the 
bright feature at 11-13 km s _1 appear on the near side of the 
peak with respect to V* (Fig. 7) at epoch 1996a, but are located 
on the far side at epochs just a few weeks later. The peak veloc- 
ity with respect to V* increases by ^ 1 km s _1 over this period. 
IK Tau also shows hints of such behaviour in Fig. 6. This spatial 
re-broadening cannot be explained by current maser theory. 

5.2. Matter-bounded masers 

The other two AGB stars, U Her and U Ori, and some IK 
Tau features, show a very loose relationship between s v and /. 
Inspection of their velocity profiles (Figs. 4 and 5) shows that the 
components in some bright features at some epochs have the in- 
verse relationship between observed size and intensity expected 
for spherical clouds. The components in other features, however, 
show very little variation with intensity, and some brighter fea- 
tures exhibit larger components than do fainter features at the 
same epoch. EHM92 explain how , in 'matter-bounded' beam- 
ing, the measured maser size is expected to be independent 
of frequency, i.e. location in the line profile. This is also con- 
sistent with the larger values of Q est obtained for U Ori than 
for any other source, although these have large uncertainties 
(Section 3.2.4). 

Matter-bounded masers emanate from regions where the am- 
plification path along the line of sight is longer than one or both 
physical dimensions of the emitting region in the plane of the 
sky. EHM92 develop a model for maser amplification through a 
flattened cloud perpendicular to the shock direction using disc- 
like masers (resembling a coin viewed edge-on), and filamentary 
masers, which can be envisaged as narrow tubes of amplification 
pointing at the observer. They explain how the apparent size of 
such masers remains similar to the true smaller dimension of the 
emission region in the plane of the sky, with little dependence on 
intensity or position in the line profile. 

Figure 2 shows that U Her and U Ori have a higher propor- 
tion of larger components than the other AGB stars, with the 
peaks of the distribution of s v comparable to / at some epochs. 
EHM92 predicted that filamentary masers would exhibit an ir- 
regular distribution of components whilst a more linear distribu- 
tion is expected for disc masers. The distribution of maser com- 
ponents in some features around U Ori (Figs. 12 and 13) and 
U Her is indeed scattered, although some of the features which 
show an increase of component size towards the line peak do 
have apparently linear structures, possibly due to velocity gradi- 
ents. 

From EHM92 equation 2.2, a masing region can behave as 
a filamentary maser if its depth along the line of sight (21 in the 
notation of EHM92) exceeds yr times its width in the perpen- 
dicular direction (2h in the notation of EHM92). U Her and U 
Ori have r 16 (Table 3), derived by assuming spherical clouds. 
t has an inverse logarithmic dependence on \T X \, which in turn 
has a logarithmic dependence on cloud depth (as explained in 
Section 3.2.3). Hence, if the depths of some clouds along the 
line of sight is greater than the measured size /, r will be slightly 
reduced. 

Taking this into account, the filamentary condition requires 
maser clouds with a depth ~ 5 times their width. If we take I as 
the width, this implies cloud depths of 12-25 AU for U Her and 
U Ori, which is only half the maser shell depth (r - n). We do 
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see some very long component series, e.g. the feature running 
from (-20, -20) to (-10, -60) in Fig. 12 epoch 1994, about 10 
AU in extent. This may be due to the effect of a velocity gradient 
which allows us to see many components across the narrow edge 
of a disc maser. 

If clouds are initially spherical, but develop an aspect ratio 
« 5, this cannot be achieved by radiation pressure alone, as ex- 
plained at the start of this Section. The other possible cause is 
shock compression and U Her and U Ori are in fact the most 
likely candidates for strong shocks reaching the H2O maser re- 
gion. 

The H2O maser shells around U Ori and U Her have volume 
filling factors 2-3 times smaller than around RT Vir and IK Tau, 
despite having a higher mass loss rate than RT Vir, which im- 
plies that masing clouds are more disrupted around U Ori and U 
Her. Samus et al. (2008) and data from the AAVSO (American 
association of variable star observers) show that U Ori and U 
Her have regular, deep stellar pulsations with amplitudes « 5 
mag, but IK Tau is a somewhat less regular pulsator, with erratic 
amplitudes. S Per and RT Vir are semi-regular variables, with 
pulsation amplitudes of < 3 and < 1 mag, respectively. Stronger 
shocks are likely to propagate out from the stellar surfaces of U 
Ori and U Her, which could flatten maser clouds. Rudnitskij et al. 
(2000) describes signs of the impact of stellar pulsations appar- 
ent from long-term spectral monitoring of H2O masers around 
U Ori. Moreover, both U Ori and U Her have undergone OH 
maser flares likely to have originated in the inner CSE (Pataki 
& Kolena 1974; Chapman & Cohen 1985; Etoka & Le Squeren 
1997). 



6. Conclusions and summary 

We measured the properties of H2O masers around S Per, U Ori, 
U Her, IK Tau and RT Vir at multiple epochs using MERLIN. 
The emission in each 0.1 or 0.2 km s _1 velocity channel is re- 
solved into milli-arcsec scale components, which form series in 
sucessive channels generally corresponding to a single spectral 
feature with a Gaussian-like profile. Many of these features have 
systematic position-velocity gradients which allow us to esti- 
mate the physical extent of maser clouds. The number densities 
required for masing and other properties show that the clouds 
are probably discrete physical entities. The average feature peak 
brightness temperatures are 10 8 ^ 7b ^ 10 12 K. 

Most of the maser features have negative excitation temper- 
atures close to zero and modest optical depths, consistent with 
mainly unsaturated masers. The excess population of the upper 
level of the maser transition is < 1 % of the combined number 
density of both states, except for RT Vir, where it reaches ~ 10% 
at some epochs. We use this to infer maser excitation tempera- 
tures of -0.1 it % k -5 K. This leads to optical depths of 8-25 
(excluding two epochs with insufficient good data). 

We find an inverse relationship between the measurable size 
of maser components in individual velocity channels and the in- 
tensity, for S Per and for all but one epoch each of IK Tau and 
RT Vir. This is broadly consistent with the relationship predicted 
by EHM92, s v oc (In/) -0 ' 5 , for amplification-bounded masing 
from uniform, spherical clouds. The brightest H2O components 
around S Per and IK Tau shrink more steeply than predicted with 
increasing intensity, suggesting that they are partially saturated. 
In such conditions the effective observed size is reduced pro- 
portionally to the size of the unsaturated core. The relationship 
between the relative size of each component and its position in 
the line profile of amplification-bounded masers is expected to 



obey \n(s v /so) °c l/2(£u/AV D ) 2 (where sq is the size of the peak 
component in each feature and s v is the size of the component at 
velocity offset Sv). We find that \n(s v /so) has a weaker relation- 
ship with Sv. Inspection of the plots of the variation of the RT Vir 
component sizes across its line profiles shows that the brightest 
components re-broaden noticeably, which cannot be explained. 

Most U Ori and U Her features and some around IK Tau 
show contrasting behaviour more suggestive of matter-bounded 
masers. The size of maser components is almost independent of 
In /, being generally larger in proportion to the sizes of features 
(/). EHM92 predicted that observed component size is indepen- 
dent of position in the line profile for matter-bounded masing. 
Many features in U Her and U Ori are likely to meet the fila- 
mentary condition, due to disc-like masers, with an aspect ratio 
^ 5. The components of U Her and U Ori have larger observed 
sizes than for the other two AGB stars, and U Ori has the largest 
beaming angles. The random spatial distribution of components 
within some U Ori and U Her features is also consistent with 
the predictions of EHM92 for matter-bounded masers. These 
sources have the deepest and most regular pulsation periods, so 
shocks could propagate more strongly into the H2O maser shell 
and flatten the clouds. 

The two scenarios, for amplification- and matter-bounded 
masers, are compared in a cartoon in Fig. 1. The upper diagrams 
show how maser propagation through a sphere will give a very 
narrow apparent spatial size at full width half maximum (left) 
but propagation down the long axis of a disc (right) produces an 
apparent size more similar to that of the emission region. The 
star names are in order of the degree to which they present the 
type of masing behaviour. These high-resolution MERLIN ob- 
servations provide direct observational tests of the predictions 
of EHM92. Most circumstellar 22-GHz H2O masers behave as 
spherical, amplification-bounded clouds. They are probably un- 
saturated, the warm conditions and small size (relative to star- 
forming regions) providing a very efficient pump which allows 
them to reach high maser optical depths. Some of the maser fea- 
tures round the Mira variables with the deepest-amplitude pul- 
sations appear to be matter-bounded, possibly emanating from 
shock-flattened clouds. 
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Appendix A: Cloud size variations and velocity 
gradients 

The observed size of clouds within each CSE varies by a factor 
of two or more (Table 2). This Appendix investigates whether 
this can be included meaningfully in our models, followed by 
an attempt to allow for the effect of velocity gradients. The end 
results are qualitatively consistent with Section 4.1.2 but the un- 
known, probably irregular, shapes of the clouds means that we 
cannot obtain precise measurements. 

The effect of variations in size of spherical clouds (with 
purely thermal velocity dispersions) was tested by rearranging 
Equation 6 in order to plot log Sv/ 1 against log(ln /). We obtained 
slopes similar to the values of a given in Table 3, or slightly 
steeper. All S Per, IK Tau and RT Vir epochs had negative slopes, 
averages -1.3, -1.4 and -0.9, respectively; U Her and U Ori had 
large scatters with positive averages. These results were of lower 
significance than the evaluation of a described in Section 4.1.2. 
It seems that additional effects are more significant than varia- 
tions in cloud size, in producing deviations from the predicted 
value of a — -0.5. This is not surprising, since the presence of 
velocity gradients is what enables us to obtain values of / close 
to the physical size of clouds. 

We tried a different approach by starting with an expression 
predicting the variation of component size with intensity as a 
function of velocity offset from the peak of the line profile for 
each cloud. Initially we retain a purely thermal velocity disper- 
sion as we assume that all clouds are in the same saturation state 
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Fig. A.l. Variation of component size across the line profile for S Per in 
1994. The size of maser components s v relative to the size of the com- 
ponent at the line centre so as a function of offset from the line centre, 
weighted by the relative effective path length, as in Equation A.2. The 
line is a least-squares error-weighted fit to this relationship. 



in a given source. The ratio of the measured size of a maser com- 
ponent in any velocity channel, to the size of the component at 
the line peak, s v /sq, should then be a consistent function of ve- 
locity offset from the line centre. We take sq as the size of the 
component closest in velocity to the centre of the spectral profile 
fitted to each feature as described in Section 2. 

Following E92, the unsaturated absorption coefficient varies 
across the line profile of amplification-bounded masers as 
kov/kqo = e^ 6v ^ v ^ where the absorption coefficient is Kno at 
the line centre, velocity Vq, and kq v is the absorption coefficient 
at velocity offset Sv. The thermal line width is a weak function 
of distance from the star, AVd oc r~ 02 (Section 3.1.1). We esti- 
mated r by solving for the 3D structure of each H2O maser shell 
as described in Murakawa et al. (2003). 

Using these relationships in Equation 5 gives an expression 
for s v as a function of location in the line profile 




(A.l) 



However, testing this relationship on our data did not show the 
expected quadratic dependence but a much weaker exponent, be- 
tween and 1 if it was meaningful at all. 

We now attempt to extend known analytic relationships in 
order to investigate an empirical relationship describing the ob- 
served maser beaming by including the effects of velocity gradi- 
ents in Equation A. 1 . 

Figs. 8 and 15 show that many clouds possess systematic gra- 
dients of velocity with position. The features in Figs. 10 A and 
B and 16 present orderly but non-linear structure. The progres- 
sion from channel to channel of component size and location can 
be followed clearly despite twists and turns (see RT Vir inset). 
Many features in IK Tau and a few in U Her and U Ori also have 
systematic gradients although the less ordered distributions seen 
in Figs. 12 and 13 are more typical of U Ori and U Her. The to- 
tal velocity spans (AV) of features are similar in all sources and 
the ratio AV/l is about 0.1-0.2 km s 1 AIT 1 for S Per (RYC99), 
and an order of magnitude greater for the smaller AGB clouds 
(B+03). 
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The angular separations between components within a cloud, 
especially close to the peak, are very small and often cannot be 
measured accurately. In contrast, the velocity of each compo- 
nent along the line of sight is precisely known. We assume that 
the clouds are spherical, with similar velocity-position gradients 
along and perpendicular to the line of sight, allowing us to use 
I and AV to represent the maximum spatial and velocity extent 
in any direction. We approximate each cloud as an ensemble of 
overlapping sub-spheres, each comprising molecules in a chan- 
nel velocity interval and each giving rise to a measured maser 
component. We treat the features making up the multiple-peak S 
Per clouds as if they were distinct, since amplification probably 
only takes place along a path length similar to / (not the entire 
depth of the cloud, since the velocity change is too great). The 
peak component samples the full depth, /, of the feature at its 
centre, but the depth of a sub-sphere corresponding to a com- 
ponent in the line wings is limited by the physical depth /„ of a 
chord through the cloud at an offset position 61 determined by 
the local velocity gradient. We assume that 61 is proportional to 
the velocity offset from the feature peak, 6v, for that component. 
This allows us to use 61/0.51 = 6v/0.5AV. By simple trigonome- 
try, (0.5/) 2 = (0.5l v ) 2 +6l 2 , leading to IJl = y/l - 6v 2 /(0.5AV) 2 . 

Equation A. 1 is derived from expressions for amplification- 
bounded masers, so it is only applicable within features where 
s v > so, excluding much of the U Her and U Ori data from this 
analysis. We only consider features with Gaussian line profiles 
fitted with better than 3<x accuracy (A% t in Table 2), and where 
6v/Q.5AV > 1 for both wings. The number of components meet- 
ing these criteria, N/3, is given in Table 3. 

We now take l v in Equation 6 to represent the path length 
producing the amplification of a component at velocity offset 6v 
of size s v and thus replace s v /so with (s v /0.5l v l)/(so/0.5ll) in 
Equation A.l. Using the more precise velocity measurements as 
surrogates for position measurements leads to 




(A.2) 



where /3 = 2. 

Table 3 gives values of ft obtained by error-weighted least- 
squares fits to this relationship for all epochs, for features where 
AV1/2 > 3crAy 1/2 .. Results in italics are of low significance and 
suggest that this model is probably not applicable to the majority 
of features for these epochs. 

All the values of /3 are roughly 50% greater (i.e. closer to 
the predicted quadratic relationship) than the equivalent values 
of the exponent obtained by fitting to Equation A. 1 . 

The 1994 S Per data, shown in Fig. A.l, give/? « 1.8, closest 
to the expected value. The other sources with significant mea- 
surements have /3 « 0.9. Many factors could contribute to the 
lower values of /3, such as deviations from spherical symme- 
try and non-linear velocity gradients. The closer angular spacing 
seen in many cases for brighter components compared with the 
displacement of the faintest from the line centre, means that 6v 
is not a perfect surrogate for 61. 



